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Research Article 

Simulation of the Behavior of Carbon 
Nanotori during Unfolding: A Study of 
Stability and Electronic Structure 


Abstract 

In this work, we present an in silico study of the dynamic behavior of a carbon nanotorus during 
the localized breaking of the C-C bonds in its atomic network. It is shown that the unfolding of a carbon 
nanotorus is accompanied by the appearance of deformation waves running along the atomic network of 
the carbon structure. We estimate energy stability of a carbon nanotorus under deformations using the 
original method for calculating the local stresses of an atomic network. A new physical phenomenon of 
the molecular current flow in the most deformed sites of the torus atomic network is discovered. 


Abbrevations 

COMPASS: Condensed-phase Optimized Molecular Poten¬ 
tials for Atomistic Simulation Studies; DOS: Density Of State; 
SCC DFTB: Self-Consistent Charge Density Functional Tight- 
Binding; HOMO: Highest Occupied Molecular Orbital 

Introduction 

Due to geometry and topology features of the atomic 
structure, the carbon nanotori are considered as prospect 
objects for many investigations. At present, carbon 
nanostructures of the given type can be synthesized by 
the method of laser growth, ultrasound treatment, organic 
reactions, gas-phase chemical vapor deposition [1-4]. Along 
with the development of synthesis technologies there are many 
papers devoted to construct the structural models of toroidal 
carbon nanostructures. Currently, there are two large classes 
of the nanotorus structural models. The first class is formed 
by pristine nanotubes with ideal hexagonal network [5]. The 
second class is formed by the nanotubes with certain amount 
of pentagon and heptagon defects [6]. Circular geometry of 
the structure and incorporation of the defects cause unusual 
electron [7-9], mechanical and magnetic properties of nanotori 
as compared to carbon nanotubes. 

The determination of thermodynamically stable configura¬ 
tion is the important problem within the construction of the 
nanotorus structural models. At present, there is no universal 
approach to numerically estimate the stability of toroidal car¬ 
bon nanostructures. For example, Chuang et al. [10], proposed 


to estimate the stability of carbon nanotori by two criteria: 
two types of curvature (the local mean curvature and the local 
Gaussian curvature) and curvature energy distribution. At the 
same time, each of the curvature parameters should be chosen 
on the basis of torus geometry characteristics, and the extre¬ 
mal values of the curvature energy do not always coincide with 
the extremal values of the Gaussian curvature. Therefore, the 
chosen parameters of the curvature can give contradictory re¬ 
sults and, thus, may distort true data about structural model 
stability of the toroidal objects. Another criterion for estimat¬ 
ing the nanotorus stability is the calculation of the carbon 
system binding energy. In the paper [11], authors performed 
a theoretical investigation of the stability for armchair carbon 
nanotori of different forms, diameter and length on the basis 
of binding energy calculation. However, according to the au¬ 
thors of this work, the suggested criterion is very sensitive to 
such parameters of atomic framework as bond length and va¬ 
lence angle that can vary significantly depending on the form 
of considered nanotorus. Thus, in the number of cases in order 
to make a right conclusion about structural stability of a to¬ 
rus it is necessary to calculate the distribution of the average 
values of valence angle and bond length for each considered 
nanotorus. In the paper [12], authors estimated the stability 
of armchair and zigzag carbon nanotori with the Stone-Wales 
defect by the value of critical diameter when phase transition 
was observed. This transition is accompanied by the trans¬ 
forming atomic configuration of the torus from the structure 
with ideal circular form to the structure with wave-like bends. 
At the same time, in this paper such phase transition was reg¬ 
istered only qualitatively by dynamic picture of the structural 
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Results and Discussion 

Molecular-dynamic simulation of the nanotorus unfol¬ 
ding 

The object of this study was carbon nanotorus (10, 6) with 
the diameter of 20 nm and thickness of 1 nm. This torus was 
obtained as a result of defectless folding of chiral nanotube with 
corresponding geometry sizes into ring. Such structures can 
be synthesized by a high-voltage pulsed discharge in ethanol 
vapor [23]. In our numerical experiment, the chemical bonds 
were destroyed on the circle of the tube in a fixed fragment of 
the torus atomic network as shown in figure 1. Such localized 
destructions can occur during the synthesis of tori when the 
chemical bonds of the structure are destroyed under external 
factors (pressure, temperature, arc discharge and others), and 
torus is, in fact, broken. Also, these destructions may be caused 
by an action of atomic force microscope tip during probing the 
surface of synthesized carbon object. We were interested in two 
questions in the course of the study of torus dynamic behavior: 
how the atomic structure of the torus would be reconstructed 
over time of the numerical experiment and how this will change 
the energy stability of the investigated structure. 

Figure 2 shows investigated nanotorus after bond breakage 
in fixed times. It can be seen from the figure that during 
the first 30 ps three sections of different form were clearly 
distinguished in the object structure: central part in the form 
of arc, and two edge parts with the form close to straight. 
Transition zones between sections form bends that spread from 
the object edges to its center like waves with velocity of 200 
m/s. After -30 ps the deformation waves reached the structure 
center and form the bend of complex form with two almost 
straight sites on its edge. After ~6o ps the central part of the 
object was straightened, and then structure gradually took the 
form of deformed tube. After the bonds were destroyed, the 
torus straightened into the tube within 91 ps. 


Figure 1 : Atomic structure of the torus at the initial moment of localized breaking 
of chemical bonds along the circumference of the framework. 


Figure 2: Step-by-step straightening of the carbon nanotorus into the tube during 
the localized breaking of chemical bonds (A) - after 11 ps (B) - after 20 ps (C) - 
after 30 ps (D) - after 91 ps. 
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changes in the atomic network. Such estimation is not physical 
and doesn’t have enough reliability. 

Prediction of stability is very important for the investigation 
of mechanical properties of toroidal structures, in particular, its 
behavior under mechanical load. There are a large number of 
papers devoted to prediction of the nanotori strength features 
under different deformations [13-13]. For example, C. Feng and 
K. M. Liew [14], studied behavior of zigzag (9,0) and armchair 
(8,8) nanotori under mechanical stretching from both sides 
using COMPASS force field. The critical displacements for the 
toroidal structure subjected to stretching on two symmetrical 
sides of the ring were obtained. For armchair tori these values 
vary from 8.16 nm to 7.9 nm depending on the structure 
diameter, for zigzag tori - from 0.68 to 2.5 nm. It was found 
that during stretching the zigzag and armchair tori took the 
form of an elongated ellipsoid. Influence of structural defects 
on the strength properties of nanotori was investigated in 
paper [13]. The behavior of armchair nanotorus (3, 5) with 
vacancy defect under stretching was modeled within classical 
molecular dynamics method. Results of numerical experiments 
have shown that the nanotorus is destroyed with a ratio of 
displacement to torus diameter larger than 50% independent 
of the defect concentration. Besides, the loss of stability for 
toroidal structure and fracture initiation were observed in 
the area of elongated areas, not in the defect areas. So, all 
above-mentioned papers did not offer an effective tool for 
prediction of destruction area in atomic network. At the same 
time, the influence of deformation on the electronic properties 
of toroidal structures is still a relevant problem. As shown 
earlier, it is possible to tune the electronic properties of such 
carbon structures as graphene nanoribbons [16], and carbon 
nanotubes [17], using deformation. 

The current paper is devoted to solving above mentioned 
scientific problems. Its main purpose is to study the dynamic 
deformation processes occurring in carbon nanotorus during 
local destruction of chemical bonds and to predict the influence 
of these processes on the electronic properties of carbon 
nanotori. 

Materials and Methods 

In order to study the nanotori behavior under deformation 
we used classical molecular dynamics method with application 
of modified Brenner potential in order to describe interaction 
between atoms. As we noted earlier, modified Brenner potential 
was successfully applied to study deformation processes of 
graphene nanoobjects under compression and stretching [18]. 
Nose - Hoover thermostat was applied to control the temperature 
of medium. All numerical experiments were performed using 
Kvazar software [19], under 300 K with modeling time-step 
of 1 fs. In order to study the changes in electronic properties 
of carbon nanotori during unfolding we calculated density 
of state (DOS) by means of self-consistent charge density 
functional tight-binding (SCC DFTB) quantum method [20], 
using the freely available DFTB+ code [21]. Using the Mulliken 
population analysis [22] we estimated valence electron charge 
density of carbon atoms during nanotorus unfolding. 
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In order to estimate the influence of deformation waves 
spreading in atomic framework on the toroidal objects stability 
we calculated field of local stresses. These calculations were 
performed by means of original method, which applied earlier 
to study the stability of graphene nanoribbons under axial 
compression [16] and bamboo-like nanotubes [17]. By stress at 
atom we mean the difference between energy density of a torus 
and energy density of a nanotube, which form the torus. The 
following algorithm was used to calculate a local stress: 

1) Optimization of initial atomic structure of a carbon 
nanotube by minimization of the total energy using 
empirical method based on the Brenner potential; 

2) Calculation of distribution of the bulk energy density on 

nanotube’s atoms; 

3) Optimization of the structure of a nanotorus formatted 

in the result of the nanotube unfolding by minimization 
of the total energy; 

4) Calculation of the bulk energy density distribution on 
nanotorus’ atoms; 


5) Calculation of the local stress field by difference of the 
bulk energy density for a torus and a nanotube. The 
local stress for atom with number i cr was calculated 
according to formula: 


0 

w, — w, 


( 1 ) 


where w t is the bulk energy density of a nanotube which 
form the torus, W t is the bulk energy density of a torus. 


Figure 3 shows fragments of atomic framework in different 
time moments with different degrees of deformation. The 
stress distribution on torus’ atoms is shown by color for 
different fragments. The sharp bends of the atomic network 
are observed during the spread of deformation wave. They 
appear for several picoseconds. In these local areas of the 
atomic framework the stress increases up to 18 GPa on some 
atoms. This value exceeds critical value for the nanotubes 
(14 GPa [17]), that’s why some chemical bonds are broken. 
Such areas of atomic framework are highlighted by red at the 
picture. Since the deformation wave moves with large velocity, 
the structure can’t relax to minimum energy and the bonds 
are restored in the next time. We suggest that this interest 
physical phenomenon is observed within whole process of 
torus unfolding. During torus unfolding the value of maximum 
stress decreases and reaches minimum in the moment of final 
torus straightening. 


The stress peaks occurred in atomic network affect the 
energy. Figure 4 shows plot of energy change during unfolding 
a chiral torus into tube. This plot shows that the sharp jump of 
energy is observed in the moment of torus breaking. This jump 
is caused by a repulsion of opened torus edges. Further, the 
total energy of the structure decreases uniformly right up until 
the final nanotorus straightening. 



Figure 3: The stress distribution of the atomic network of a nanotorus during 
straightening at fixed instants of time (A) - 0 ps (B) - 1.5 ps (C) - 7.5 ps (D) - 15 
ps; (E) - 22.5 ps (D) - 45 ps (G) - 91 ps. 



Figure 4: Dynamics of the change in the potential energy of a nanotorus during 
straightening. 


Prediction of changes in the electronic properties of a 
chiral nanotorus during unfolding 

In order to identify the effect of the deformation wave 
on the electronic properties of a carbon torus, we calculated 
DOS. Figure 5 shows the DOS of the nanotorus at fixed time 
moments. In these moments, we observed the sharpest bends 
of the atomic network of toroidal structures. The time of 10 
ps correspondWs to the case when two symmetrical sharp 
bends were formed in the atomic network of a torus. The time 
of 22.3 ps corresponds to the case when one bend of complex 
shape was formed in the atomic network of a torus. As can be 
seen from the graph, there are no visible changes in the DOS 
distribution for the deformed nanotorus at the transition from 
two sharp bends of the framework to one bend of the complex 
shape formed in the region of the closing deformation waves 
moving towards each other. 
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The DOS distributions were calculated for the original 
defectless chiral torus and the corresponding chiral tube in 
order to estimate the changes in the electronic properties of 
a torus during the propagation of deformation waves along 
its framework. The calculation results are shown in figure 6a. 
As can be seen from this figure, the DOS distributions for the 
nanotube and torus have a pronounced similarity. However, 
in the electronic structure of the torus there is an important 
change in comparison with the nanotube structure, namely, 
the appearance of an energy gap of 0.5 eV. Figure 6b shows 
an isolated fragment of the DOS distribution near the Highest 
Occupied Molecular Orbital (HOMO) level, clearly illustrating 
the appearance of an energy gap in the electronic structure of a 
defectless chiral nanotorus. Also, this figure shows a fragment 
of the DOS distribution of a straightening nanotorus at the time 
of the largest bending of the atomic network. Analyzing the 
presented results of the DOS distributions, we can conclude that 
as soon as the torus loses the integrity of its atomic framework 
and ceases to be a closed structure, it loses semiconductor 
properties and becomes a conductor. 

Further, we estimated valence electron charge density 
of carbon atoms during nanotorus straightening using 
the Mulliken population analysis [23]. In the course of the 
study, it was found that the charge transfer between the 
carbon atoms of the deformed framework was observed with 
the nanotorus straightening. The maximum values of the 
electron charge transfer (~ o.ose ~ o.09e) were found in the 
regions of the largest bends of the atomic structure of the 
torus. The maximum stresses are also concentrated in these 
structural regions. This fact is due to the previously established 
relationship between the mechanical characteristic and the 
electron density of the atomic framework [16]. According to 
our previous results obtained for carbon nanostructures, the 
rehybridization electronic orbitals will be stronger in the areas 
of the largest curvature of the surface. That is why the largest 
charge transfer is fixed in the region of the largest curvature of 
the atomic grid of the torus as a consequence of deformation. 
Using the maximum value of the electron charge flow, we 
numerically determined the amount of current flowing in the 



Figure 5: The DOS distribution of a nanotorus at fixed times, in which the sharpest 
bends of the atomic network were observed. 



Figure 6: The DOS distribution for a nanotube and a nanotorus formed from it: a) a 
general view; b) a distribution fragment illustrating the change in the DOS near the 
HOMO level (~ 6.2 eV) for a nanotube, torus formed from it, and straightening torus 
at the time of maximum inflection. 



Figure 7: Dynamics of the current change in a carbon nanotorus during its 
straightening. 


carbon nanotorus during its straightening. Figure 7 shows the 
dynamics of the current change during relaxation time of the 
toroidal structure. 

It can be seen from the above plot, the maximum value of 
the current (~ 1.1 nA) is reached 7.3 ps after localized breaking 
the chemical bonds on a circle of the nanotorus. At this time, 
the atomic structure of the torus is characterized by the 
presence of four symmetrical bends of the framework. As the 
torus is straightened, the current begins to decrease, reaching 
a minimum at the instant of the disappearance of the sharp 
bends of the atomic network. 

Conclusion 

The phenomenon of formation and propagation of 
deformation waves in the toroidal nanostructures was revealed 
during the study of dynamic behavior of a chiral nanotorus 
at localized destruction of the chemical bonds of the atomic 
network. These waves appear after breaking torus and lead to 
its deformations and numerous local bond breaks. Since the 
deformation wave moves along the atomic framework of the 
torus with a high velocity (~ 200 m/s), the toroidal structure 
does not have time to relax to a minimum of energy, and the 
chemical bonds are restored at the next time. In this regard, 
the atomic framework of the torus remains defect-free after 
being unfolded into the tube. The results of calculation of the 
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local stress distribution indicate that the straightened torus 
is characterized by high energy stability. The time of the full 
unfolding process for a chiral torus with a diameter of 20 nm 
is ~ 91 ps. 

Analysis of the calculation results of the DOS distribution 
showed that at the instant of localized breaking of chemical 
bonds along the circumference of the tubular framework, 
the torus loses its semiconductor properties and become 
a clearly expressed conductor. During a study of electron 
charge distribution at different instants of the unfolding of 
a nanotorus, we found a new physical effect, consisting in 
that the molecular current flows in the torus atomic network 
sections that have undergone the greatest deformations. 
The maximum value of the current depends on the degree of 
deformation of the atomic framework. 

Thus, the obtained results allow us to consider carbon 
nanotori as a promising material for nanoelectronics that is 
resistant to mechanical stresses and deformations. 
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